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ABSTRACT 

In many recent papers concerned with providing an expla- 
nation for the geomagnetic anomaly, favorable agreement with 
measured data has been obtained from the equations of motion 
for electrons and ions when used with an empirical boundary 
condition, whereas poor agreement has resulted from attempts 
to numerically integrate the diffusion equation derived from 
the equations of motion. We have been able to demonstrate that 
the deviation of the commonly employed diffusion equation is 
based on assumptions concerning the equations of motion which 
are inconsistent with the observed distribution unless an 
additional constraint equation is also applied. Unfortunately, 
application of the constraint equation also provides incorrect 
results. We find, therefore, that none of the currently used 
forms of the diffusion equation takin g geomagnetic control into 
account provide correct results because they are based oh improper 
physical models for the ionosphere. 

Since the equations of motion do provide a favorable des- 
cription for the geomagnetic anomaly, we have studied the pos- 
sible physical models leading to the form of the equations used, 

! 

and found that although field aligned diffusive equilibrium 
provides the correct form, a more reasonable assumption 


I 


concerning electron and ion collisions with neutrals also leads 
to the same result. We have also been able to provide a more 
realistic theoretical description of the geomagnetic anomaly by 
employing an analytic form for the boundary condition which is 
more accurate with measurement than those previously used. 

Finally, by combining the equations of motion for neutrals, 
electrons and ions, we have been able to predict geomagnetic 
control for the neutral atmosphere in the lower F region of the 
ionosphere, although the exact shape of this distribution’ is 
unknown. 
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GEOMAGNETIC CONTROL OF DIFFUSION IN THE UPPER ATMOSPHERE 


Introduction 

In recent months it has become increasingly evident that 
a certain amount of confusion exists in the understanding of 
thejhasic physical mechanisms governing diffusion and the 
existence of the geomagnetic anomaly in the ionosphoreT^ This 
apparent confusion arises by comparison of the work of Chandra 
(1964), (to be referred to as C-I) in which it is shown that 
the assumption of ambipolar diffusion along a field line can- 
not lead to geomagnetic control of the charged particles, and 
such papers as Goldberg and Schmerling (1962, 1963), (to be 
referred to as GS-I and GS-II) and Goldberg, Kendall, and 
Schmerling (1964), (to be referred to as GKS) in which this 
process does appear to produce geomagnetic control of the 

charged p&rticXe density in the ionosphere. It is our contention 

'' ... . 

that these two apparently opposite view points do not contradict 
each other and that the confusion is almost entirely one of term- 
inology usage and comprehension of the fundamental physics gov- 
erning the derived equations. 

The purpose of this paper is to describe and resolve the 
confusion which exists in thefield at this moment, and then to 
point out the new problem with which we must contend in order 
to derive and apply the diffusion equation to ionospheric prob- 
lems correctly. In addition, a section will be devoted to an 
improved theoretical .description of the geomagnetic anomaly, by 
using an analytic expression for the vertical electron density 

distribution at the equator which is more in accordance with 

measurement than the simple Chapman type distrubution employed 
in GKS. **> 
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Fundamental Equations and Definitions 

The major cause of confusion appears to lie in the appli- 
cation of two phrases, viz. ambipolar diffusion and diffusive 
equilibrium . Let us investigate and discuss each of these terms 
to determine how loose usage of them has led to the current 
problems of misunderstanding. 

In the normal sense, ambipolar diffusion refers to a plasma 
in which the negative (electrons) and positive (ions) charges do 
not move independently due to the influence of the electric field 
caused by their Coulomb interactions. In this medium, the elec- 
trons and ions drift in pairs and this motion of electron- ion 
pairs is referred to as ambipolar diffusion. The condition for 
ambipolar diffusion in a neutral plasma is thereby 

v 0 - V. = V Cl) 

e 1 

where v is macroscopic velocity and the subscripts e and i refer 
to electrons and ions respectively. When ; 

v - 0 (2) 

the condition for diffusive equilibrium is then satisfied. 

f 

The implications of (1) are quite straight-forward, as 
shown in C-I. ‘In an isothermal atmosphere and in the presence 
of magnetic field, this requires ?X(vxB) “0. In particular, 
the assumption of field aligned plasma diffusion (vxB » 0) can 
only be satisfied for a trivial case, v = 0, resulting in hydro- 
static distribution of electron density independent of geomag- 
netic latitude. 
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On the other hand, favorable comparison between Alouette 
topside sounder measurements and theoretical calculations of 
the geomagnetic anomaly has been obtained in GKS by assuming 
conditions of ambipolar diffusion and diffusive equilibrium 
along field lines, thereby indicating a possible conflict with 
the results in C-I. The problem resolves itself once one investi- 
gates the meaning of ambipolar diffusion and diffusive equili- 
brium in the GKS sense. 

Let us first write the general equations of motion for 
neutrals, electrons and ions, respectively, where the subscript 
n refers to neutrals. Following C-I: 


n in nr . a ^ xi • m « m ^ 

rfi % + ttt V <v-i> - -% + vv 

m e + v ®i + ®n 


n.-m m. 

■ i e i 


Wn 


vei (v e -vj) + • e -- e ven (v fe 

»+ m i m e + % 
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n.ra.in 
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vei (Vi - v e ) + 


n i m i m n 
m i + 


vin (v t - v n ) 


-vp^ + n^tug + en^ (E + v^ x B) 

where n is number density, is the collision frequency be- 
tween the and particle, m is mass, p is pressure, g is 
gravitational acceleration, e is the absolute value of electron 
charge, E is electric field, and B is magnetic field. In 

- - 2 - - 
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writing equations (3} - (5) it is assumed that a . 

n k n Z 

In the following we assume that the plasma is in a quasi - 
neutral state 


n ft# n. *» N 
e i 


and;, the electrons, ions and neutrals obey the ideal gas law in 
the ionosphere, 

p -n.kT. (7) 

J J 3 

when k is Boltzmann's constant and T is temperature. Further- 
more, we assume thermal equilibrium, i.e. 


T “ T “ T « 

e i 

Then, 


p e - - p 


In addition, we assume for simplicity that . . 


v n «0 


(10) 


Then summation of (4) and (5) provides 


m n m e N V ® n in i N 

ven v e + 

m e + mn m ± + m^ 


(U) 


vin - -2?p + N (m e + g + tJ x B 
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where xJ is current density^ defined as 


xJ - Ne (v. - v ) 
i , e 


(12) 
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Since we are investigating ambipolar diffusion and diffusive 
equilibrium in the GKS sense, it is desirable to write this 
equation in component form along a field line as 


m n m e _ 

veo v 

m + in ^ e 
e n 
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B 
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( 13 ) 
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when B is a unit vector in the direction of the magnetic field. 
Let us write (13) in more familiar form by using 

m e << “n U 4 ) 


and defining the scale height of the ionizable constituent as 


where 


H. 


kT 


i m.g 

Also, for convenience, we make the approximation 


(15) 


^ ~ m i 


(16) 


Then 


m e v en -* 
-V- v e 


A Uli 

B + T v in 


_ A 

v.- . B 


-w <- *s+ JL>- b 




Finally, we write 


(17) 


A * A 

B - -(i r sinl + i cos I) 


(18) 


A • * 

where i r and ig are unit vectors in the r and 0 directions and 
I is the magnetic dip angle, reckoned positive when the north 
seeking pole of the needle points downward.. Now, if we treat 
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ambipolar diffusion in the GKS sense, we simply imply that the 
electron and ion velocity components in the field direction are 


equal, i.e, 


«♦ ^ -» A 

v . B » v. . B - v 

e 1 x i 


(19) 


Applying (18) and (19) in (17), we obtain 

v ,i - ^ Uin i d iS+ij + I M 3 

* N 3r Nr 3e 


( 20 ) 


where 


H v - m e en 
2 


m i v in 


( 21 ) 


Assuming that 


mV « m. V 
e en x m 


( 22 ) 


we may write 


jmv » m i v in 


(23) 


because of (14) . 


Equation (20) is a familiar result derived in such papers as 
Kendall (1962) and GS-II. However, it is clearly not the result 
of ambipolar diffusion, which is given by (1), but instead, the 
result of a statement concerning the field line components of 
electron and ion velocities given by (19). 


If we now demand 


( 24 ) 
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which is the statement implying diffusive equilibrium along a 
field line in the GKS sense, we obtain the familiar equation 

sin I (IM+ 1 ) 4 . cos I _ a (25) 

n ar + 7 + ~W VS 0 

which can also be written as 
1 dN 1 - 

N ctr + -207- 0 (26) 

provided we recognize that r and 0 are not independent in (25) 
but related by the dipole field condition 

r » r Q sin^ 0 (27) 

~ t>. 

(an^tan 2 m 2 co/ ^ (28) 

It is evident that (26) can only be treated in this total deri- 
vative form if the integration is carried out along the field 
line. 

Statements concerning the components of vectors in a par- 
ticular direction, rsuch as (19), do not imply any conditions 
on the total vector. As a result, (25) has not required the 
assumption of any restrictions on the behavior of the velocity 
components normal to the field lines. 

Equation (26) has been the basis of describing geomagnetic 
control in the upper F-region in GKS paper. Although the results 
of this paper appear justified on the grounds that diffusive equi- 
librium occurs along a field line, it is nevertheless undesir- 
able to apply a concept whi^ch is of purely hypothetical nature. 

We now investigate other assumptions to find a more realistic 

/ 

justification for (25). 
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Consider equation (17). If we apply (19), we obtain 


m e v en 


2 ^^ + ,T i v in v i - -IcT (— + %—) . 
2 e xl ^4 ln N 2 h/ 


(29) 


We find two ways in which the right hand side will approach zero. 
The first imposes a new ...ondition on the velocities, viz: 


V e 7 ~ m i V in v i 

. 2m e v en 


v = -m f v. v,. 
e 1 . a- 

11 2m v 11 
e en 


iZ --- 






(30) 


a result which, although possible, appears rather unlikely since 

it would require a very special condition that the electron 
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velocity be of the order of 10 times greater and in the opposite 
direction than the ion velocity. 

However, if we can demand that the collision frequencies 
between electron and neutral and ion and neutral be sufficiently 
small so that the drag forces arising due to collision be neglig- 
ible as compared to the' pressure gradient, gravity and Lorentz 
forces, it is possible to derive equation C3l 6) without imposing 
any restriction on the velocities of electron and ion. We be- 
lieve that this assumption is more realistic in the upper F-region 
where the gyro- frequencies of electron and ion are much greater 
than their corresponding collision frequencies. 

Although the collision frequency assumption is physically 

I 

more desirable, it prevents us from obtaining a simple expres- 
sion for'Wjjj^ or ^ . Instead, we must return to the original 


U ' 
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equations of motion, (4) and (5), and solve for 4^ and ^ expli- 

n/ 

citly, as has been carried out in the appendix in C-1. Unfor- 
tunately this introduces a very serious complication in the work 
because of the difficulty of eliminating electric field from the 
expressions of v e and v^ without making specific assumptions 
about the relationship between v e and v^. The implications of 
these assumptions will be discussed in the latter part of the' 
paper. In the following section we proceed to discuss the physi 
cal implications of equation (26). 


The Electron Density Distribution with the Effect of a Variable 
Scale Height 

Equation (26) can be integrated along, a field to provide 
the general solution 


'r cot^ 0 


N(r, 0) f(r , tt/ 2 ) eV 2 H. 


(31) 


However, if we treat T and constant but recognize that g is 
proportional to 1/r^, H. is proportional to r^, and we obtain 


NCr.O) - f(r ot 


U» p 
2 M t - (A ) 


(32) 


In both cases, f(r , tt/ 9 ) is an arbitrary function of height 
at the equator which cannot be determined by the equations of 
motion from which (31) or (32) are derived. The function 
f(r o ,n/ 2 ) must therefore be given as a boundary condition in 
this problem and can only be determined empirically or by use 
of additional equations governing the physics of the problem. 
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Since (31) or (32) depend exclusively upon the equations 
of motion, it appears that; an additional equation such as the 
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continuity equation, should lead to. the desired boundary con- 

■ . ia\ i«' v ‘ K: i. 

— : dition. Unfortunately, as we -have-di3cussed-earlier , the deri- 

vation and solution of the continuity equation depend upon a 
knowledge of E. Thus, the complexity of the problem becomes 
quite formidable and it is difficult to anticipate a simple 
method of solution at this time. 

Instead we depend upon an empirical type boundary condi- 
tion, which may very well be the solution of the correct conti- 
equation, to derive the explicit form of the electron density 
distribution. 

The incorporation of a Chapman like boundary condition 
in (31) leads to the results obtained in GKS. Since such a 
boundary condition can only be considered as a rough approxi- 
mation to the shape of the actual vertical electron density 
distribution at the equator, it is desirable to employ an ana- 
lytic boundary condition which more closely resembles the true 
height profiles. Chandra (1962) has proposed a modified form 
of the Chapman function which includes the effect of variable 
scale height and which is found to fit the measured vertical 
-«• distribution for electron density at mid -latitudes far more 
accurately than the simple Chapman form. We assume here that 
such a function also describes the electron density distribution 
at the equator. Thus we may write 


f(r 0 ,rr/2) - N r exp \ “ r j- — r~° 

m 9 ^J.1 "• a 6X11 \~ a ^ r ° 


r o " r mo 


- exp 


2H 0 .1 - a exp a ( r o 

2H c 

/ ? 


' r mo) 
l o 4 


( 33 ) 
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where H~ is the scale height of A ionic constituent and N- is 
o o a . r mo 

the value of electron density at" equatorial height r mo . The 
parameter which is a measure of departure from simple Chapman 


function^is defined as 


a 


- H(r^) 
o mo 


(34) 


where H(r mo ) is that value of H(r rt ) at r. 


•mo* 


Also, r n is 


understood to be the radial height specifically at 6 - ~/2. 
If (32) is substituted into (31), we obtain 
^ — r csc^ 0 - 


N(r, 6) 
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r exp H. 
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(35) 


Equation (35) then provides a general expression for the electron 
density at all heights and co -latitudes provided that we are 

in a region where the effect of collision can be neglected. 

. - -A 

The variation of/N r with co-latitude are shown in 

Ant \ _ 

,i.x Wm. a ^ 

(• u- •• *- s v- 

figures 1-3 for different values ° a * ^mi and H 0 . ^ For— tne 
purpose— of— -numer ical com putation- the radius of;\earth is taken 
to be 6370 km, and it is assumed that H 0 • H^. It is seen that 
the basic features of geomagnetic anomaly are unaltered by 

I 

changing the various parameters. - Further, in equation (35), if 

we identify the term H 0 Lfl-a exp — (co-see^ 0 - r mo )"l with- — 

w t 1 ' 2H ft Ojc. J R 


^ / 
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in;,GKS paper, it is seen that since:, a is positive kH^ > l/ l for all 

heights This explains why Kh> 1 gives better fit with the experi- 
mental data in the -latter paper. 

Problems Involved in the Derivation of the 33iffusion Equation 
In the previous sections we have- seen how the equations 
of motion for electrons and ions are sufficient to obtain a theo- 
retical description of the electron density distribution in the 
topside equatorial region of the ionosphere under equinox con- . 
ditions. This has required us to make certain assumptions con- 
cerning collision frequencies or velocity components along field 
lines and also forced the application of an empirical boundary 
condition at the equator. In order to produce the empirical 
boundary condition theoretically and also obtain a solution which 
is valid in both the topside and bottomside equatorial F region, 
it is necessary to turn to the continuity equation for additional 
information. Using the explicit expressions for velocity which 
are derivable from the equations of motion, it is then possible 
to derive the diffusion equations associated with the ionosphere. 

If we simply require total ambipolar diffusion to occur in 
the ionosphere so that v is independent of the electric field 
explicitly, and furthermore demand that in all regions concerned, 
the motion along field lines are much larger than the drifts 
normal to field lines, we must . then invoke the additional con- 
straint . equation 

- A 
v x B “ o 

iM. 


(36) 


This leads to a hydrostatic distribution of electron density 
which doesn't agree with measured results, as has been demon- 
strated in C-I. A second approach (Kendall, 1962, and GS II) 
is the assumption that ambipolar diffusion exist only along field 
lines (see equation (19)). Thus, if we assume that the parallel 
velocity components of electron o t ion velocity are equal and 
much greater than either of the unequal normal velocity components, 
we can write (20) as a good approximation for the entire velocity. 


In mathematical notation, we have 




(37) 


where and v^ are the perpendicular components of electron 
and ion macroscopic drift; velocity respectively. This implies 


that 


Vn « V 


(38) 


and provides us with a velocity expression independent of elec- 
tric field. The general contention here has been that although 
(37) implies (38) so that we can write 


Q - L = v . N v « V .Nv xi 


(39) 


as the steady state continuity equation, where Q and L are pro- 
duction and loss respectively, we no longer need invoke the 
equation of constraint since (36) is not exactly true. We can 
then substitute (20) into (39) and obtain the well known form 
of the two dimensional diffusion equation with no additional 
equation of constraint. 

We wish to discuss this approach by first questioning the 
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validity of (37), and then demonstrating that even if it were 
true, C28) cannot imply (39) without the additional inclusion 
of (36). This should demonstrate that the field line ambipolar 
diffusion approach with neglect of the normal velocity components 
is completely identical to the total ambipolar case in which 
the major component of drift velocities are assumed to lie along 
field lines. Thus, the results of the two approaches will be 
identical, leading to the conclusion that ambipolar diffusion 
in which the major component of drift velocity lies along a field 
line cannot be the correct physical model to describe the equa- 
torial electron density distribution in the F region of the 
ionosphere . 

Let us first consider (37). We have already seen that 
the assumption of diffusive equilibrium along a field line 
( = o) leads to a correct description of the electron density 
in at least the topside region of the ionosphere. If this is . 
the true model of the physical situation, then it is inconsis- 
tent with (37) and we cannot expect any results obtained using 
(37) to provide us with correct results concerning this region. 
If, on the other hand, the neglect of momentum transfer terms 
instead of diffusive equilibrium^ can be attributed to small 
collision frequencies^ (37) need not be violated. This might 
be a further justification for validating the collision frequency 
assumption instead of the diffusive equilibrium model. Unfor- 
tunately, as we approach the equator, we see from (20) that v^ 
approaches zero since both sin I and 3N/50 approach zero. The 
latter condition is based strictly on the empirical condition 


of symmetry about the equator. We therefore find that no matter 

how sm all v and v may be, there will always be a region about 

e ± 

the equator in which (37) does not apply unless 

v = v. = o (40) 

e x X x 

which is identical, to the equation of constraint, (36). 

We now return to the second question. That is, even if (37) were 

true, which could still be possible provided v ^ v. , can we 

6 il x \\ 

describe the electron density distributions in the entire region of the 
ionosphere by (39) ^ We note that 

V . Nv = v • ■+ N V . v = (v^ Hrv^) • + K 7 . (v^-v v (l ) (4l) 

since 

v = v 4< v, (4fi) 

it X 

Now, in order to write (39)* we must demand that 

\y . VH + H V . | > > |v L . VK -j- N V . vj (43) 

We first note that the geomagnetic anomaly is symmetric about the equator 
during equinox. This implies that » 

l 1 1 ) 
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11m v u * VN = lim v B . = 0 (4-4- ) 

6 - rr /2 9 - tt/2 

so that in some region about the equator, (43) becomes 

|n V . V I > > |w V . + v x . ?n| (45) 

Although (45) could be true .for certain special cases, there is no 
a priori guarantee that (45) will be implied by ( 37 ) in general without 
the additional condition that Vj_ = 0 . Thus, if we are to write (39) as 

a direct and general implication of ( 37 )* we are once again forced to 

— ♦ A 

employ the constraint equation v x B = 0. 

We cannot state that (40) holds in a very small region about the 
equator so that its effect outside this region can be neglected. The 
geomagnetic anomaly itself is a second order effect and we cannot 
expect to reproduce it by neglecting second order terms which are res- 
ponsible for its existence. 

We therefore find that if ambipolar diffusion exists in the 

ionosphere, and if it is restricted to the field line direction, we 

cannot assume ( 37 ) without imposing an additional constraint equation. 

Furthermore, (37) does not generally imply ( 39 ) in the ionosphere with 

or without ambipolar diffusion unless the constraint equation is also 

employed. However, since (37)* (59)* and the assumption of ambipolar 

diffusion along a field line do not provide 

1 
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ttt 7 the correct description of equatorial electron density, <N, 

| * f 

j we must conclude that these assumptions are not valid in a theory 

i i 

j leading to a description of the electron density distribution : 

I ! 

in the ionosphere. f 

Kendall (1962), and Rishbeth, Lyon and Peart (1963), have | 

| i 

I attempted to numerically integrate (30) derived from (20) and 

! (37), without invoking (26). They have been unable to obtain 

the correct description of the geomagnetic anomaly and have 

therefore concluded that diffusion may not be a very important 

i; 

physical process governing the empirical distribution of elec- 
tron density. However, on the basis of the discussion presented 
in this section, it now appears that the physical assumptions j 

used in deriving the form of (39) used in their work may not 
be valid, which simply implies that the diffusion equation is 

far more complicated than originally believed. Since (19) and [ 

f 

(37) are no longer valid, we can no longer equate electron and i 

f- ‘ j 

ion velocities to eliminate electric field. Insted, we must j 

■j 

write separate continuity equations for electrons and ions and j 

i: 

describe the behavior of electric field before it is possible [ 

to obtain the correct theoretical description of the geomag- | 

netic anomaly from the diffusion equation. j 

It may appear that the results presented in GS II are j 

also not valid for the reasons discussed above. However, a [ 

I, 

closer inspection of GS II shows that no new information was 

| 

obtained from the solution of continuity equation than that \ 

already available from the equations of motion. The equation [ 

i 1 

\S ' I 


discussed in GS II was simply 

' V- N 5-J-- 0 

?Nv , « 0 (46) 

, 11 

I 

0 

Where the explicit production and loss terms were neglected in 
obtaining the series solution. Furthermore, as shown in GKS, 
the equation of motion leading to (26), whether derived assuming 
v ii“ ° or n^hing assumptions concerning collision terms^has the 

identical solution to that obtained from (46) in GS II. (For D 

*(? * 
<v xl obviously cannot give any new information.) This 

explains why the empirical boundry condition was necessary to 

obtain a non-arbitrary solution using (.39) in GS II. 

— ‘ We should point out however, that solutions of (4'6) 

making use of explicit production and loss terms would not give 

correct results at^L'e^st' in the equatorial region for the reasons 

discussed in this section. 


The Distribution of the Neutral Atmosphere 


In this section we will show that when the. drag forces 
are not negligible, as would be the case in the lower F-region 
and E-region, it is possible to study the behavior of the neutral 
atmosphere without imposing any restriction on the velocities 
of the various constituents. To obtain the necessary starting 
equation,' we first sum (3), (4), and (5): 


P'i 

- 


- v (p e + Pi + p n ) + (n n m n + N.) g + Zf x B “ 0 


(47) 


where we have once again u£ed (14). The component of (47) along 
the direction of magnetic field is then 
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C" 7 Cp e + + P n ) + (m n n n + m i N) 13 • B b 0 


(48) 


Comparison of (47) and (48) shows that the net force due to 
pressure gradient and gravity of all particles is perpendicular 
to the magnetic field and balanced by a current flow force. Next, 
using (7), (8), (9), (18) j and (28), we have 


S-*(n n + 2N) tan 8 ^ (n n 4- 2N) n n N 

+ + 4 - - 


(49) 


H n Hi 


where ^ is the scale height of the neutral atmosphere, 
Since N << n n and Hi « H n , we can write 

n N ru. N ' r! ru+ 2N 

— + — w — ^ 4 — 

Hn Hi H n H,, Hn 


(50) 


whereASuffix on n has been dropped for simplicity, 
becomes, in total derivative form, 


Then, (49) 


d(n 4 - 2n) + n 4 - 2N 


(51) 


Integration of (51) can be obtained along the field line giving 


-J'r £ 


n(r, 0) « n + 2N « g(r 0 , Tr/2) e o 


(52) 


where g(r Q , tt/2) is an arbitrary function of height at the equa- 
tor and r Q is defined in (27). If we now demand that the radial 
distribution of the neutrals obey the normal hydrostatic law at 


the equator, so that 


g(r Q , tt/ 2) - n 00 e r 


r r o 


(53) 



where n QO is the neutral number density at height r OQ on the 
equator, then 


-r dr . 

" n oo e r °° ^ » 


(54) 


a result which is entirely independent of 9. If, on the other 
hand, g(r Q , tt/ 2) is perturbed in any manner from the exact hydro- 
static equilibrium case, we will obtain a distribution for n 
which does depend on 9. The origin of this angular dependence 
on the neutrals may seem somewhat surprising until we realize 
that in selecting a functional form for g(r Q , n /2), any devia- 
tion in the equatorial neutral distribution from hydrostatic 
equilibrium must arise due to collisions between neutrals and 
geomagnet ically controlled charged particles. Thus, if the 
collisions between neutrals and charged particles are sufficiently 
large to make the momentum transfer forces between charged and 
neutral particles important, the neutrals will begin to tend 
toward the angular distribution of the geomagnet ically controlled 
particles. This can, also be seen from (3), where it is obvious 
that we will not obtain the exact hydrostatic distribution in 
a region when the terms on the left hand side become important. 

On this basis, we might expect to observe angular variations 
of the neutral distribution in the bottomside regions of the 

ionosphere where charged -neutral particle interactions become 

/ 

important. 

t 

Conclusions 

From the discussion and results of this -paper, we have 
, shown the following: 


i 


1. From the equations of motion, it is possible to 
derive an expression for electron-density distribution along 

a field line either by assuming diffusive equilibrium along the 
direction of the magnetic field or by neglecting the drag forces 
arising due to collision. The latter assumption appears to be 
more realistic in the topside of the ionosphere. In either case, 
it is necessary to assume a distribution at a certain point in 
the radial direction in order to obtain a complete description 
of the electron density distribution. 

2. We have provided a more accurate formula for the repre- 
sentation of the equinox geomagnetic anomaly than that produced 
in GKS. Since the empirical boundary condition equation^.has 
been shown by Chandra (1962) to fit nearly all vertical profiles 
of electron density measured to date, we can safely assume that 
the proper selection of parameters in this formula will lead 

to a reasonable reproduction of the anomaly in any region of 
the ionosphere where either simultaneous diffusive equilibrium 
of ions and electrons occurs, or where interactions of neutrals 
with charged particles are small because of infrequent collisions. 

3. The derivation of the diffusion equation commonly 
employed to obtain a theoretical description of the electron 
density distribution in the F region ionosphere depends upon 
the assumptions v e * v,* » v^ or v,- . Attempts to numerically 

1 1 ii 

integrate the diffusion equation derived in this manner have not 
described the geomagnetic anqmaly with the correct behavior. 

We have now been able to show that the assumptions made in the 
deviation are inconsistent with experimental data unless the 



‘additional equation 'of constraint v x B =* 0 is also employed. 

Since the numerical integration results are based on an 
incomplete derivation of the diffusion equation, it is under- 
standable why such results are in poor agreement with measurement. 
Furthermore, we find that the inclusion of the constraint equa- 
tion leads to a horizontally stratified electron,^ distribution 
which also f is incapable of describing the geomagnetic anomaly. 

We must therefor, conclude that the diffusion equation derived 
in the above manner is based on assumptions .which do not fit 
the true physical description of the ionosphere and^lead to the 
neglect of the very physical parameters responsible for the geo- 
magnetic anomaly. 

4. A study of the neutral atmosphere distribution has 
led us to the conclusion that geomagnetic control of neutrals 
occurs in any region of the ionosphere where interactions of 
neutrals with charged particles become important, since this 
is most likely to occur in the lower F region of the ionosphere, 
we suggest that such geomagnetic control of the neutrals might 
be observable in this region. 
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